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Abstract We developed a self-consistent conserving pseudo particle approxima- 
tion for the Anderson impurity model with finite Coulomb interaction, 
derivable from a Luttinger Ward functional. It contains an infinite series 
of skeleton diagrams built out of fully renormalized Green's functions. 
The choice of diagrams is motivated by the Schrieffer Wolff transforma- 
tion which shows that singly and doubly occupied states should appear 
in all bare diagrams symmetrically. Our numerical results for Tk are in 
excellent agreement with the exact values known from the Bethe ansatz 
solution. The low energy physics of non-Fermi liquid Anderson impurity 
systems is correctly described while the present approximation fails to 
describe Fermi liquid systems, since some important coherent spin flip 
and charge transfer processes are not yet included. It is believed that 
CTMA (Conserving T-matrix approximation) diagrams will recover also 
Fermi liquid behavior for Anderson models with finite Coulomb inter- 
action as they do for infinite Coulomb interaction. 

1. INTRODUCTION 

Anderson impurity models have been of considerable interest recently 
for studying the transport in mesoscopic structures [1] as well as for de- 
scribing strong correlations on the lattice in the limit of large dimensions 
[2]. In both cases the doubly occupied impurity site, which is character- 
ized by an additional energy cost [/, is essential for the correct physical 
description. 

In terms of pseudofermion (/o-) and slave boson operators (a, h) the 
Anderson model is defined by the Hamiltonian 
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+ E ^kic^h^fa + <^ctj+^a + h.c.) + i7a+a; 

where is the energy of the local d level, U is the Coulomb interaction 
and are the hybridization matrix elements. A physical electron in 
the local level is created with the operator da = b'^ fa + af^^a, where 
b represents the empty level and a the doubly occupied one (in the 
following we will call a and b the heavy and the light boson, respectively). 
For the representation to be exact, the local operator constraint Q = 
J2a fa fa+b^b+a^a = 1 must be fulfilled at all times. This is guaranteed 
by a suitable projection technique [4]. 

2. METHOD 

It is well known that the Non-Crossing Approximation (NCA) gives 
the correct Kondo temperature Tk in the limit of large Coulomb repul- 
sion U while it fails to recover the correct Kondo temperature for finite 
U. The reason may be seen from the Schriffer- Wolff transformation [3] 
where the effective coupling constant has two equally important terms 
for U ~ -2€d > 

J=-— + — -7 (1.1) 

where just one of them survives in the large U limit. In Eq. (1.1) the first 
(second) term is generated by exchange of a bare light (heavy) boson. 
This implies that the light (5) and heavy (a) bosons should appear sym- 
metrically in all bare diagrams. In order to generate the corresponding 
diagrams of the Luttinger-Ward functional one may start with the NCA 
diagram for the U ^ oo case (first diagram in Fig. la)) and replace the 
light boson line by a heavy boson, exchanging the conduction electron 
lines appropriately. The result is the second diagram. This is not suffi- 
cient, however, because the diagrams should be symmetric in a, b on the 
level of bare diagrams. One may show that this exchange generates an 
infinite series of diagrams, each containing a circle of pseudoparticle lines 
with one light boson and n heavy bosons (or vice versa) dressed by con- 
duction electron lines spanning at most three pseudoparticle lines (see 
Fig. la)). The corresponding self-energies are obtained by functional 
differentiation and are shown in Fig. lb) in terms of vertex functions 
(the dash- and cross-shaded half-circle areas). 

The Green's function of the physical d electrons can be calculated from 
the exact relation that can be derived with the help of the equation of 
motion 

= hk'GU^) + Eg.-,GL-(-)^^G,(a;)FiG;-, -,(-). (1.2) 
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Figure 1 a)The Luttinger Ward functional, b) the Bethe-Salpeter equations for the 
T- matrix vertices (dark and light shaded areas), c) the Contributions to the irreducible 
self energies for pseudo-particles as well as for conduction electrons are shown. 

The local conduction electron Green's function can be expressed in terms 
of the local conduction electron self energy Sc (see Fig. Ic)), 

Grf(^) = (1.3) 



3. NUMERICAL RESULTS 

The Kondo temperature was extracted from the width of the Abrikosov- 
Suhl peak in the local electron Green's function. The values we get from 
the present approximation are in excellent agreement with the exact val- 
ues known from the Bethe ansatz solution, as can be seen in Fig. 2. The 
simple Non-Crossing Approximation (the first two terms in $ in Fig. la) 
gives Kondo temperatures several orders of magnitude too low. Includ- 
ing one additional diagram does not increase Tk sufficiently {UNCA) 
(here we disagree with the corresponding claim in Ref. [5]). An infinite 
series of skeleton ladder diagrams is needed to get the correct scaling of 
Tk with the Coulomb energy U and the local electron energy e^. 
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Figure 2 Left: Kondo temperature for various parameters (Bd and f/). Solid lines are 
values of Tk from the Bethe-ansatz solutions. Right; Local electron spectral function 
calculated with NCA, UNCA and the ladder approximation. The width of the peak 
is proportional to the Kondo temperature. 

4. SUMMARY 

In summary, we have shown that an infinite number of skeleton dia- 
grams is necessary to recover the correct dynamic energy scale [Tk) in 
the Anderson impurity model with finite Coulomb interaction energy. 
We found a systematic way to obtain those diagrams and have shown 
that the approximation is derivable a from Luttinger-Ward functional 
and therefore conserving. Numerical results show excellent agreement 
with the exact values of Tk- 
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